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Abstract. - The persistence length of semiflexible polymers and one-dimensional fluid mem- 
branes is obtained from the renormalization of their bending rigidity. The renormalized bending 
rigidity is calculated using an exact real-space functional renormalization group transforma- 
tion based on a mapping to the one-dimensional Heisenberg model. The renormalized bending 
rigidity vanishes exponentially at large length scales and its asymptotic behaviour is used to 
define the persistence length. For semiflexible polymers, our result agrees with the persistence 
length obtained using the asymptotic behaviour of tangent correlation functions. Our defini- 
tion differs from the one commonly used for fiuid membranes, which is based on a perturbative 
renormalization of the bending rigidity. 



Introduction. ~ Thermal fluctuations of two-dimensional fluid membranes and one- 
dimensional semiflexible polymers or filaments are governed by their bending energy and can 
be characterized using the concept of a persistence length Lp. In the absence of thermal fluc- 
tuations at zero temperature, fiuid membranes are planar and filaments are straight because 
of their bending rigidity. Sufficiently large and thermally fluctuating membranes or filaments 
lose their planar or straight conformation. Only subsystems of size L Lp appear rigid and 
maintain an average planar or straight conformation with a preferred normal or tangent direc- 
tion, respectively. Larger membranes or fllaments of sizes L ^ Lp, on the other hand, appear 
flexible. In the "semiflexible" regime for which L is smaller than or comparable with Lp, the 
statistical mechanics is governed by the competition of the thermal energy T and the bending 
rigidity k. Experimental values for the persistence length of one-dimensional biological fila- 
ments vary from 50 nm for double-stranded DNA [T], 10 /xm for actin filaments up to the 
mm-range for microtubules [2_. The persistence lengths of two-dimensional fluid membranes 
composed of lipid bilayers are typically much larger than experimental length scales. 

For semiflexible polymers with one internal dimension, Lp is usually defined by the char- 
acteristic length scale for the exponential decay of the two-point correlation function between 
unit tangent vectors t along the polymer. A continuous model for an inextensible semiflexible 
polymer of contour length L is the worm-like chain (WLC) model [4] . In the WLC model the 
polymer is parametrized by its arc length s (0 < s < L) and the polymer contour is completely 
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determined by the field t(s) of unit tangent vectors. The Hamiltonian is given by the bending 
energy 

^{t(s)} [ ds {dstf, with t^{s) = 1. (1) 

^ Jo 

where kq is the (unrenormalized) bending rigidity of the model. For a WLC embedded in d 
spatial dimensions, the tangent correlation function is found to be [S] 

(t(s)•t(s')) = e-l^-^'l/^^ withL,^^^. (2) 

where T is the temperature in energy units. 

For fluid membranes with two internal dimensions, the analogous quantity is the correlation 
function of normal vectors. An approximate result has been given in Ref. [5], but a rigorous 
treatment is missing due to the more involved differential geometry. Surfaces cannot be fully 
determined by specifying an arbitrary set of normal vectors, but have to fulfill additional 
compatibility conditions in terms of the metric and curvature tensors, the equations of Gauss, 
Mainardi and Codazzi, which ensure their continuity. Implementations of these constraints 
lead to a considerably more complicated field theory than ([1} describing a two-dimensional 
fluid membrane in terms of its normal and tangent vector fields [7] . 

For fiuid membranes, an alternative definition of the persistence length Lp has been given, 
which is linked to the effect of K-renormalization. The mode coupling between thermal shape 
fluctuations of different wave lengths modifies the large scale bending behavior, which can 
be described by an effective or renormalized bending rigidity k. The renormalized k has 
been calculated using different perturbative renormalization group (RG) approaches [5Hl3). 
The results are still controversial: Several authors [SHU] find a thermal softening of the 
membrane with increasing length scales, but differing prefactors, whereas Pinnow and Helfrich 
[13] obtained the opposite result. Furthermore, different definitions of the persistence length 
are considered in these approaches: In Refs. [5HTT]. Lp is identified with the length scale, 
where the renormalized bending rigidity k vanishes, while Helfrich and Pinnow defined Lp via 
the averaged absorbed area [51[TB]. 

In this Letter, we obtain an exact real-space RG scheme for the bending rigidity of a 
semiflexible polymer or a one-dimensional fluid membrane, which allows us to define the 
persistence length as the characteristic decay length of the renormalized bending rigidity. 

In principle, a perturbative result for the effective k can be deduced from the RG analysis 
of the one-dimensional nonlinear ct- model, which is equivalent to the WLC Hamiltonian ([Ij. 
After a Wilson- type momentum-shell RG analysis, one obtains the effective rigidity (see, e.g., 
Refs. [n]) 

T d- 2 f 1 It . 9 , 9 J 
kq tt a Aq J J 

which depends on the momentum A. The parameter kq — k(Ao) is the 'bare' coupling taken 
at the high momentum cut-off Aq = 7r/6o, which is given by a 'lattice spacing' or bond length 
bo- Using also A = n/iwe obtain the renormalized k = k(^) as a function of the length scale £. 
Following the procedure previously used for membranes, the persistence length can be defined 
via 

K{Lp) = and, thus, Lp ~ J^Y' 

For the case of the polymer in the plane, the Hamiltonian simplifies to a free or Gaussian field 
theory such that n — kq is unrenormalized to all orders in kq/T and, thus, Lp as defined via 
K{Lp) = would become infinitely large. 



k{A) _ Ko 
T ~ Y 
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A similar perturbative momentum-shell RG procedure is possible in the so-called Monge 
parametrization of a weakly bent semiflexible polymer, analogous to the RG analysis for two- 
dimensional membranes [3]. Then the polymer is parametrized by its projected length x 
with < a: < L^:, where is the fixed projected length of the semiflexible polymer while 
its contour length becomes a fluctuating quantity. The renormalized n — k,{IIx) becomes a 
function of the projected length scale ix, which complicates a comparison with the result 
([3]) , which was derived in an ensemble of fixed contour length. Using the analogous criterion 
K,{Lp) = we obtain Lp ~ 27r^Ko/(3(i — 1)T within the Monge parametrization. 

A comparison of the RG results from the non-linear cr- model, see eq. and the one 
obtained in the Monge parametrization with eq. ^ shows that the RG results for the persis- 
tence length Lp are not compatible with the definition using the tangent correlation function. 
This raises the general question which of the definitions should be preferred. 

In this work we concentrate on a discrete description for semifiexible polymers, which is 
equivalent to the one-dimensional classical Heisenberg model. The advantage of this model is, 
that the K-renormalization as well as the tangent correlation function are exactly computable 
in arbitrary dimensions d. Consequently a direct comparison of the persistence length deter- 
mined via K-renormalization and via the tangent correlation function is possible. We introduce 
this model in the next section. The K-renormalization is carried out in a similar fashion as 
is commonly used for Ising-like spin systems. In contrast to the nonlinear cr-model, we find 
nontrivial results for K{t} both in two and in three dimensions. As expected for an exact 
result, is always positive and approaches zero only asymptotically. Finally, we analyze 
the large scale behavior of k.{1) leading to a power series of exponentials with the same decay 
length as obtained for the tangent-tangent correlations. We define this length scale to be the 
persistence length of the polymer. 

Theoretical model. - A discretization of the WLC Hamiltonian ([1} should preserve its 
local inextensibility. In addition, we want to use a discretized Hamiltonian which is locally 
invariant with respect to full rotations of single tangents - in addition to the global rotational 
symmetry of the polymer as a whole. A suitable discrete model is an inextensible semifiexible 
chain model as given by [T5] 

^{t^}-^^*' (l-t.-t,_i), witht^^l, (5) 

with M bonds or chain segments of fixed length 69. The semiflexible chain model is equivalent 
to the one-dimensional classical Heisenberg model (except for the first term, which represents 
a constant energy term) describing a one-dimensional chain of classical spins. 
The partition sum reads 



Zm = 



n / dt, exp[-H{t,}/r] =\X{ dtA \{_^ r,,,_i , (6) 



where we have introduced the transfer matrix 



Ti,i^i = exp [-^0(1 - U ■ ti-i)] , with Kq = Ko/boT. (7) 

We can parametrize the scalar product of unit tangent vectors using the azimuthal angle 
difference A9i^i-i as tj • tj_i = cos(A(?i.i_i). Then the transfer matrix can be expanded as 

00 

Tm-i= E Xili\Ko) = e-^oi„,{Ko) (8a) 



m— — 00 
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in two dimensions and 

oo 

T,.,_i = ^(2Z + 1)a[°V,(cos A^,,,_i), A[°'(i^o) - 



2A7 



°Il + l/2iKo) 



(8b) 



in three dimensions, where Ik{x) denotes the modified Bessel function of the first kind and 
Pi{x) the Legendre poiynomials [TB]. In the following, the sums J2m=~oo ioT d — 2 and 
^^(,(2/ + 1) for d = 3 are abbreviated by J^^''- 

For simplicity, we restricted our analysis to d = 2 and d ~ 3 spatial dimensions, but 
our results can easily be generalized to arbitrary dimensions d: The transfer matrix is then 
expanded in Gegenbauer polynomials and the eigenvalues aJ"'' are proportional to modified 
Bessel functions Ii^ci/2-i{Ko)- 

The partition sum and tangent-tangent correlations may be calculated exactly for open 
and periodic boundary conditions as was done, e.g., in d = 3 by Fisher [T7] and Joyce jl8j . 
For arbitrary dimension d, the tangent-tangent correlation for open boundary conditions is 
simply 

... . - . L/bo 

(9) 



which reduces in the continuum limit of small bo or large Kq to 

Renormalization Procedure. - The real-space functional RG analysis for the semiflexible 
chain ([5|) proceeds in close analogy to the one-dimensional Heisenberg model |19| and similarly 
to the Ising-like case where the t^'s are confined to discrete values [5D]. Similar real-space 
functional RG methods have also been used to study wetting transitions or the unbinding 
transitions of strings [HK^^]. In each RG step, every second tangent degree of freedom is 
eliminated. We introduce a general transfer matrix 



Ti,i^i = exp [h[ti ■ ti^i,K)] 



(10) 



where h = h{u, K) defines an arbitrary interaction function depending on the scalar product 
of adjacent tangents u — ti ■ ti_i and the parameter K. We start the RG procedure with an 
initial value K = Kq and an initial interaction function h{u, K) = —K{1 — u), see eq. ([7]). 
Also for an arbitrary interaction function h{u, K) we can expand the transfer matrix in the 
same sets of functions as in ([S]), which defines eigenvalues X„i — \m{K) in 2d and A; = Xi{K) 



in 3d. Initially, these eigenvalues are given by X„i{K) = \m\K) and Xm{K) = X\"'{K), see 
©. 

Integration over one intermediate tangent t' between t and t" defines a recursion formula 
resulting in a new interaction function h' = h' {u, K) and an energy shift g' by 

exp[/i'(t • t", K) + g'{K)] ^ j dt' exp[ h{t ■ t', K) + h{t' ■ t", K)], (11) 

where the energy shift g' is determined by the condition that h'{l,K) — h{l,K) = 0, i.e., the 
energy is shifted in such a way that the interaction term is zero for a straight polymer. This 
leads to 

(12) 



exp[g'{K)] = J dt exp[2 h{t ■ t' , K)]. 



The recursions ([TT]) and ([T2l) are exact and can be used to obtain an exact RG relation for 



the eigenvalues A^^"* after N RG recursions. 



A 



(N+l) 



AN) 



Si: K"' 



(13) 
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In general, the new and old mteractions h'{u^ K) and K) will differ in their functional 
structure. Thus the renormalization of the parameter K cannot be carried out in an exact and 
simple manner as for one-dimensional Ising-like models with discrete spin orientation |20) . The 
only fixed point function of the recursion (jlip is independent of u, i.e., h*{u,K) = because 
of h*{l,K) = 0. This result, together with the condition h'{l,K) — 0, which is imposed at 
every RG step, suggests that the function h'{u, K) can be approximated by a linear function 

h'{u, K) ~ -K'{K) (1 - u) for w = t ■ t" ~ 1, (14) 

as long as the scalar product w = t • t" is close to one, i.e., sufficiently close to the straight 
configuration. This approximation should improve when the whole function h! {u, K) becomes 
small upon approaching the fixed point h*{u,K) — after many iterations, i.e., on large 
length scales. Using the approximation (fT4|) . K'{K) is defined by the slope of h'{u,K) at 
u = 1, 



d 



K'(K) = V' ^ = —eyii,\h'(u,K)] , with u = t • t". (15) 

du u=i du '■ ■' 11=1 

Equivalently, one could expand the explicit expression for h'(x, K) given by (|lip and the right 
hand side of (fTl|) for small tangent angles and compare the coefficients. In order to extract the 
renormalized bending rigidity k' from the result for K' , one has to take into account that K' 
also contains the new bond length b' = 2b, which increases by a factor of 2 at each decimation 
step. Therefore, 

k'{K) ^2bTK'{K). (16) 

Using this procedure we can calculate the renormalized bending rigidity kn after N RG re- 
cursions in 2d and 3d starting from the exact RG recursions (|13|) for the eigenvalues. Inserting 
the renormalized eigenvalues into taking the derivative according to (fT5| and applying 
the rescaling (ITB)) finally yields the result 



^{i:l''[A<»)(A-.)f"A«}/{5:r[A™(A.)f"}, (IT) 



Kq Kq 



with = and A^^"^ = ^n{n + 1). In the following we will interpret kat as a continuous 
function k{£) of the length scale £ by replacing the rescaling factor 2^ = bN/bo by the 
continuous parameter i/bo. 

Persistence length. - The sums in the expressions for the effective bending rigidity (jl7p 
can be computed numerically. Fig. [T] displays the results for k{£)/ko as a function oi £/ba for 
Kq = 1000 and in 2d and 3d. The value Kq = 1000 is appropriate for a semiflexible polymer 
with kq/T = 10 /im and a bond length fog = 10 nm, which is close to experimental values for 
F-actin [2l[3]. For DNA, appropriate values are kq/T ~ 50nm and bo ~ 0.3nm and, thus, 
Ko ~ 150. 

As long as £ is small, k decays almost linearly in d = 3, which is also in qualitative 
agreement with the result from the RG of the nonlinear u-model. For d ~ 2 the decay is 
much slower at small length scales, but, in contrast to the non-linear cr-model where k is not 
renormalized. This qualitative difference is due to the following important difference between 
the Heisenberg and the nonlinear a-model: Parametrizing the WLC model ([Ij via tangent 
angle leaves only quadratic terms cx (Adi^i-i)'^ , whereas the discrete semiflexible chain ([S]) 
gives terms cx 1 — cos(A0i^i_i), which represent the full expansion of the cosine and obey the 
local invariance under full rotations. 

As £ increases k{£) approaches zero only asymptotically. Therefore, the definition of the 
persistence length as length scale where the renormalized k vanishes, K{Lp) = - which is 
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Fig. 1 - Renormalized bending rigidity k{£)/ko as a function of the length scale £/bo — 2^ for 
Ko ~ 1000 for d = 2 (o) and d = 3 (•) according to the recursion relation (|17|l . The lines show the 
asymptotic behavior for kq/T and ^ <C kq/T according to eqs. ((TSJ and (PU)) . respectively. 



usually used for fluid membranes - would always give an infinite result. We propose not to 
ask at which length scale the renormalized k reaches zero, but rather how it reaches zero. For 
I > boKo = kq/T the sums in (IT71) converge fast and one has to include only the first few terms 
for accurate results. In fact, one can replace the Bessel functions contained in the eigenvalues 
A^°\ see dH), by their asymptotic form Ik{x) « (x/27r)"^/^ exp[a; — (fc^ — l/4)(2x)~^] for large x 
|16j . This is justified for sufficiently large Kq > 100, which is fulfilled by semiflexible polymers 
like F-actin {Kq ~ 1000) or DNA {Kq ~ 150). Using this asymptotics we find {xI^\Kq)Y/''« ~ 
^-m^i/2boKo for d = 2 and {\f\Ko)Y/''° - e-Ki+i)e/2boKo ^j. = 3. Moreover, we may 
expand (IT71) as a power series in e~^'^^'^° and obtain 



for d = 2, 



k{£)/ko « {£T/ko) {ie-'^/^" ~ 9e-^'T/-° + 426-^^^/''° - . . .) for d 



= 3 



(18) 



The characteristic length scales in the expansions are 2kq/T in d = 2 and kq/T in d — 3, which 
are, therefore, a natural definition for the persistence length Lp. For general dimensionality 
d, the exponent of the first term is determined by the order of the Bessel function appearing 
in the eigenvalue. Thus the RG calculation leads to a persistence length 



2ko 



T{d-lY 



(19) 



which agrees exactly with the result ^ based on the tangent correlation function. 

Our definition based on the large-scale asymptotics of the exact RG fiow is qualitatively 
different from the definition dH used in perturbative RG calculations. While the result ^ 
from the nonlinear cr-model is only valid for small length scales t <C kq/T, where k{£) is close 
to Ko, the expansion (IT^ describes the region £ 3> kq/T. Indeed, taking the expansion of ([TT]) 
for £ < Ko/T, that is 



k{£) /nn^l- (87r2«;o/^T)e-2'''«o/^T ^ o{£e-'^''''''>/^'^) for d = 2, 
k{£)/ko w 1 - {£T/6ko) - 0{£^T^/kI) for d = 3 , 



(20) 
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and defining Lp by tlie exponential decay length, in two dimensions, respectively, by the linear 
term in three dimensions leads to a persistence length, which is considerably bigger than the 
value found above. Hence it is not surprising that the two computations and yield 
different results. The slow exponential decay in the expansion ()20p for d = 2 is reminiscent of 
the non-renormalization of n in the non-linear cr- model, see eq. ([3]), and leads to a "plateau" 
in the numerical result for k{£)/ko for £ <C kq/T in Fig.[T] 

Conclusion. - In conclusion we have presented a definition of the persistence length Lp of 
a semiflexible polymer or one-dimensional fluid membrane based on the large scale behaviour of 
the RG flow of the bending rigidity k, as obtained from a functional real-space RG calculation. 
Our result (|19|) for Lp generalizes the conventional definition based on the exponential decay 
of a particular two-point tangent correlation function and gives identical results for Lp, thus 
justifying past experimental and theoretical work based on this conventional definition. The 
RG flows P?)l or allow us to follow the behaviour of a semiflexible polymer from a 
stiff polymer on short length scales to an effectively flexible polymer on large length scales 
quantitatively as a function of the length scale. On large length scales, our functional RG 
gives qualitatively different results from perturbative RG techniques, which have been used 
for the closely related problem of fluid membranes [5HT5] and which we also applied to the 
one-dimensional semiflexible polymer. The generalization of our renormalization approach to 
two-dimensional fluid membranes is complicated by the more involved differential geometry of 
these two-dimensional objects and remains an open issue for future investigation. In analogy 
with our results for one-dimensional polymers and fllaments as described here, we expect 
that a proper deflnition of the persistence length for two-dimensional fluid membranes again 
requires the aymptotic RG flow on large length scales which remains to be determined. 
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